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£ : ABSTRACT 

in 

CN ■ The paper is devoted to the analysis of particle acceleration in Gamma- Ray 

^ ■ Bursts (GRBs) and its radiative consequences. Therefore we get on one hand 

constraints on the physics and on the other hand possible signatures of particle 
acceleration that could be recorded by the new gamma ray instruments. In a 
recent paper we have shown that UHECRs can be generated in GRBs even with 
conservative assumptions on the magnetic field and the scattering capability of 
its perturbations, provided that a suitable relativistic Fermi process is at work 
during the so-called " internal shock" phase. In this paper we extend the analysis 
of the consequences of these assumptions to the whole prompt emission of both 
electrons and protons. Indeed, assuming that the magnetic field decays in 1/r 2 
and that the scattering time of particles is longer than the Bohm's assumption, 



in particular with a rule derived from Kolmogorov scaling, we show that the five 
following events naturally happen with no other parameter adaptation that the 
intensity of the magnetic field, that turns out to be subequipartition: i) UHECRs 
can be generated with a sufficient flux (~ 1 km _2 yr~ 1 ) within the GZK-sphere 
to account for the CR-spectrum at the ankle (in the previous paper, we showed 
that the associated p7-neutrino emission is tiny), ii) The peak energy of the 
gamma spectrum around 100 keV, namely the so-called E pea k, is conveniently 
explained, iii) A thermal component below the E pea k is often unavoidable, iv) 
The cosmic rays could radiate gamma rays around 67 MeV (in the co-moving 
frame, which implies ~ 20 GeV for the observer) due to 7r°-decay and a low 
energy neutrino emission (around 0.2 GeV) associated to neutron decay and also 
neutrinos of energy between 5 and 150 GeV from muon decay (as predicted in 
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the previous paper), v) The UHECRs radiate high energy gamma rays between a 
few hundreds of MeV and 10 GeV (taking the pair creation process into account) 
due to their synchrotron emission with a sufficient flux to be observable. 

Subject headings: gamma-rays:bursts - particle acceleration - neutrinos 

1. Introduction 

Whereas the afterglow emission of Gamma-Ray Bursts (hereafter GRBs) has been well 
investigated and provided a strong support to the "fireball model" (Rees & Meszaros 1992), 
the origin of the prompt emission of is not well established yet. The most widely accepted 
scenario explaining the burst emission is the internal shock model (Rees & Meszaros 1994): 
according to this model, the prompt gamma-ray emission results from either the synchrotron 
emission of accelerated electrons or the inverse Compton scattering off the synchrotron pho- 
tons, the relativistic electrons being produced by the collision of relativistic shells. But, 
even if the prompt emission spectrum is correctly fitted by the empirical Band Function 
(Band et al. 1993), no physical process satisfactorily explains it especially below the peak 
energy (E pea k). More recently, the possibility of a thermal component as the low energy 
part (< 100 keV) of the prompt emission spectrum had been studied (Ghirlanda et al. 2003) 
and successfully compared to observations. In previous papers (Gialis & Pelletier (2003), 
Gialis & Pelletier (2004)), we were interested by the issue of the particle acceleration in the 
internal shocks and by the resulting cosmic ray generation: our results significantly changed 
the usual interpretation of the Fermi acceleration (Waxman 1995) and this led to propose 
an additionnal Fermi process to reach the goal of the UHECR production. In this paper, we 
intend to emphasize the consequences of the particle acceleration on high energy gamma-ray 
emission and we give an interpretation of the energy spectrum of the prompt emission. We 
also predict both the production of lower energy photons in the UV-X-ray range and a very 
high energy emission (between a few hundreds of MeV and 10 GeV, free from opacity effect) 
which should be detectable, for instance, with the GLAST experiment. 
The paper is organized as follows: in the section 2, we briefly describe the fireball dynamics 
from the central object to the deceleration radius. We also calculate the photospheric radius 
and we determine some radiative parameters we need for this paper. The section 3 summa- 
rizes our previous results (Gialis & Pelletier 2003) concerning the cosmic ray acceleration in 
GRBs and extends the study to the electron population. We analyse the consequences of the 
particle acceleration on high energy emission in the section 4. The last section of this paper 
is our conclusion on the new results that we have obtained and which could be confirmed by 
forthcoming experiments such as HESS2 and GLAST. 
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2. Dynamics and radiative parameters 



2.1. Outline of the fireball dynamics 



In this subsection, we briefly summarize all the results we need for this paper that de- 
scribe the fireball dynamics (Meszaros et al. 1993) before the deceleration radius where an 
external shock starts. 

We choose to describe the outflow with some simplified hypothesis: the wind flow can 
be considered to be a set of discrete shells which are successively emitted with an energy 
E s = E/N s , where N s is the total number of shells and E the total energy released by the 
fireball (Daigne & Mochkovitch (1998)). We assume that the total energy radiated in gamma- 
rays E 1 is a sizeable fraction of E. According to the observations, this flow is collimated 
with an average opening solid angle, Q, of about 47r/500 (Frail et al. (2001)). The wind flow 
duration, namely At w , provides with an interval of shell number which is 1 < N s < cAt w /r 
where r ~ 10 7 cm is the size of the central object. In a primeval stage of a shell expansion, 
the radiative pressure gives the temperature, T, of the completely optical thick plasma which 
is mainly composed by electron-positron pairs and by electrons and protons beyond the pair 
annihilation radius. This temperature can initially be defined by 



where E s is the energy of an emitted shell and a s = 7.56 x 10~ 15 (c.g.s) is the Stefan constant. 
For an energy E s of the order of 10 51 erg, the temperature at ro is about a few MeV. 
A characteristic baryon loading parameter, 77, is defined as the ratio between E and the 
baryon rest mass energy : 



where is the total baryonic mass ejected. The value of the baryon loading parameter is 
usually considered between 10 2 and 10 3 , in order to solve the "compactness problem". 
In a first stage, the ejected shell follows an adiabatic expansion in the surrounding medium 
and its internal energy is progressively converted into kinetic energy. Hereafter, we define 
the stationary frame as the rest frame of the central object. In the stationary frame, we 
can define a radius, r s , where the kinetic energy of baryonic matter reaches its saturation 
value. At this moment, the Lorentz factor T of a given shell is close to 77 which constitutes 
an average value. This last point is important to understand the internal shock model. 
In the same frame, the shell thickness, Ar, remains constant and equal to r until the 
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broadening radius r b > r s (Goodman (1986), Meszaros et al. (1993)). Beyond this radius, 
the shell thickness becomes Ar ~ r/2T 2 ~ r/rj 2 , so that the broadening radius is about rf r$. 
In the co-moving frame of a shell, entropy and energy conservation before r s give the evolution 
of the Lorentz factor T, the temperature T and the shell thickness, namely AR = T Ar: we 
have T oc r, T oc r _1 and AR = r/r^. Thus, the saturation radius r s is equal to rjro- 
Beyond r s , the Lorentz factor T remains constant and the previous evolution laws become 
r ~ rj and T oc r~ 2 / 3 . We can also write 



At last, the shell thickness in the co-moving frame is such that AR = rjro before r b and 
AR = r/r] beyond. 

The internal shock model has been designed (Meszaros et al. 1993) in order to account for 
rapid variability observed in GRB light curves and which can reach the millisecond. The 
internal shock model scheme is the following one: let us consider two shells leaving the 
central engine separated by a time interval At, respectively with the Lorentz factors F 1 and 
T 2 such that r 2 > IV A collision occurs at the date: 



Assuming an instantaneous shock pulse, some time spreading, At s = t c /2T 2 , is observed. 
Thus, the shortest variabilities which are observed (~ 1 ms) will be such that tb/2T 2 ~ r /2c. 
These correspond to typical time scale associated to the size of a black hole of a few tens 
solar masses (namely r /c). Such a first collision takes place around the distance r b . Longer 
variations correspond to collisions at a more remote distance until a maximum distance 
determined by the duration of the flow At w . This maximum distance is r max ~ r b c At w /ro, 
with cAt w /r ~ 3 x 10 3 (At„,/ls), which gives a proper length of the flow in the co-moving 
frame £ = f3 cT At w . The duration of the flow during the internal shock phase is therefore 
At max ~ ( r b/ r o) At w ~ rf At w . This phenomena is observed during a time interval shortened 
by the propagation effect, namely At ofcs = (1 — /3) At max ~ At max /2rf ~ At w . The previous 
value of r max is not far from the deceleration radius, r d) of the shells which is about 10 16 cm. 
Also, the Fermi acceleration of particles, which is usually considered (Waxman 1995), takes 
place in the range that extends from r b to r d , namely the internal shock phase. 




2/3 



keV 
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2.2. Radiative parameters 



We have seen that, in a primeval stage, the ejected plasma is optically thick with respect 
to the Compton scattering. Using the results of the previous subsection, we propose here to 
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determine the photospheric radius and we will define some radiative parameters . 

First, it can easily be checked that a typical shell width AR becomes smaller than the 
flow transverse radius after a short while, when r > rj^ir / '4f2r which is comparable to r s . 
It will turn out that the photosphere is located at a much larger distance for large enough 77 
and therefore the opacity of a shell is determined by its width. Assuming the temperature 
is such that 7 e hv <C m e c 2 where 7 e is the average electron Lorentz factor (fy e ~ 1 beyond 
the pair annihilation radius), the optical depth can be defined by r* = <jTn e AR with <jt is 
the Thomson cross section. Because of the plasma neutrality, we have n e ~ n p , so that the 
co-moving electron density can be written 

n ^ sE 

r Vtr 2 ARm p c 2 ' 

where £ s = 1/N S . Thus, the optical depth is 



(5) 



T * = TT~2 — ■ 6 

\ir z m p c A r\ 

We can define a critical value for rj such that the photospheric radius (t* = 1) is located at 
r b , where shock acceleration starts. This critical value r)± is given by 

-™(^) ,/5 (^)"fe) 1/5 - 

For a GRB with N s ~ 100, ^ ~ 1100 and drops around 450 for a long GRB with about 10 4 
shells. 

Thus, we can express the photospheric radius r* as 

'n \ 5/2 

V -) ■ (8) 

For usual values of 77 and according to Eq. (7), we conclude that the photospheric radius 
can be over r b , so that the internal shocks start accelerating particles in an optically thick 
plasma. We will analyse some consequences in the next sections. Also, beyond r*, one can 
consider that photons and electrons decouple and, if r > r&, electrons can be accelerated via 
the Fermi acceleration in the internal shocks. 



Considering the resulting black body emission at the photospheric radius, Eq. (3) gives the 
temperature which is such that : 

T , =0 . 37 (im(j-)-^(^y W3 keV . (9 ) 

\5MeV J V300/ \rj J K ' 
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This result must be compared to the following one: we have x L y = rf 2 Qrl<rT* where x is 
defined as the ratio between the average black body luminosity, namely L^, and the GRB 
gamma-ray luminosity L 1 < E/At w . We deduce the expression 

The comparison between Eqs. (9) and (10) leads to a ratio L^jy easily reaching a few 
percent: 

h± = 3.8 x io- (im' (w*) ( y y (^y y ~ 5/3 . ud 

y \5MeVJ \2 x 10- 3 / \10 51 erg. s- l ) V300/ \rj J K J 

Thus, the result is a thermal component in the GRB spectrum which can be observed, before 
an higher energy emission, in the range 10 keV — 180 keV in the observer frame as we see 
in Fig. 1. Some authors have already considered this possibility (see e.g Ghirlanda et al. 
(2003)) which seems to be consistent with observations. Moreover, we note that Eqs. (9) 
and (11) indicate a more important thermal component around 100 keV (with rj = 400) for 
a high number of shells because of a small i]* parameter. 



3. Particle acceleration in the internal shock phase 

3.1. Proton acceleration and cosmic ray generation 

The Fermi acceleration (first or second order) in the internal shock model is usually 
considered (Waxman 1995) as mildly or sub-relativistic with a characteristic time propor- 
tional to the Larmor time (Bohm scaling). However, in a previous paper (Gialis & Pelletier 
2003), we have shown that this assumption is not realistic regarding the magnetic energy 
depletion time. Moreover, the Fermi acceleration time depends on the mean free path, £, 
of the particle in an irregular magnetic field. This length depends on two other lengths, 
namely the Larmor radius, ri, and the correlation length, £ c : for a turbulence spectrum of 
magnetic perturbations in a power law of index (3, the following law, which is known in weak 
turbulence theory, has been extended in the regime of strong turbulence and large rigidities 
such that r L < i c (Casse et al. 2001): 

where r/ t = > • The Bohm scaling I ~ r L , which holds for electrostatic turbulence, 
does not apply with purely magnetic irregularities on large scale; no theory nor numerical 
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Fig. 1. — Top: Variation of the black body radiation temperature, measured in the observer 
frame (T b. s — V^*) as a function of the baryon loading parameter 77. Bottom: Variation of 
the ratio L bb /L^. Left: the GRBs with N s = 100 and 77* = 1100. Right: the long GRBs with 
N s = 10 4 and 77* = 450. 

simulation has confirmed Bohm's conjecture. The Bohm estimate corresponds only to the 
specific case where the magnetic field is totally disorganized and the Larmor radius as large 
as the correlation length which is not the case in GRBs. The realistic situation is the 
following one: assuming that the correlation length is of the order of a shell thickness, 
namely AR, the characteristic acceleration time is t acc — nti where ti is the Larmor time 
and k oc (r^/Ai?) 1- ^. According to a Kolmogorov scaling with f3 = 5/3, and defining kq as 
the ratio of the acceleration time over the Larmor time for a Larmor radius that equals the 
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correlation length of the magnetic field, we have: 

where e is the energy of a proton and assuming that the magnetic field strength decreases like 
r~ a . Although unproved, this scaling is more reasonable and this conservative assumption 
will lead to sensible results. 

Comparing this time with the expansion time, t exp = r/cr], we have shown (Gialis & Pelletier 
2003) that GRBs are unable to produce UHECRs with this acceleration process because of 
a strong expansion limitation in energy beyond r b which is 



e^l.SxlO'^j ^J^^-j GeV, (14) 

This limitation, measured in the co-moving frame, is more severe than the synchrotron one 
and suggests that we have to consider another type of process to achieve high energy. In 
a recent paper (Gialis & Pelletier 2004), we have proposed a Fermi acceleration process 
resulting from scattering off relativistic hydromagnetic fronts at the very beginning of the 
internal shock phase. The efficiency of this process, as shown by numerical simulation, 
is sufficient to generate a sizeable fraction of UHE cosmic rays. Also, we found that this 
scenario could constitute a very interesting additionnal acceleration process which stretch the 
cosmic ray distribution tail obtained with the usual Fermi acceleration process. Moreover, 
we showed that, for a magnetic field strength decreasing like r -2 , the high energy cosmic ray 
population is such that: 

— oc 7 , (15) 

where N* is the number of UHE cosmic rays and 7 their Lorentz factor, this spectrum 
extending over 4 decades from 10 7 to 10 11 GeV in the observer frame. 

Assuming such a spectrum for the whole proton population in the co-moving Lorentz factor 
range [1, 10 9 ], we will have: 

N pl ~ 2 , (16) 



d 7 

where N p is the total number of protons released by GRB. For N p ~ 10 51 , the number of 
UHECRs (7 > 10 8 ) generated by GRB is about 10 43 . Considering a GRB rate of about 
1 per 10 6 Mpc~ 3 yr _1 (see e.g. van Putten & Regimbau (2003) or Frail et al. (2001)), we 
deduce that, in the GZK-sphere of ~ 1 Gpc 3 , GRBs release 10 46 UHECRs yr" 1 . Because of 
the intergalactic magnetic field, this UHECR population is almost isotropized, so that we 
can observe 10 46 /47r (1 Gpc) 2 ~ 1 UHECR km _2 yr _1 which is about the observed flux at the 
ankle in the UHECR spectrum. For a better estimate of this flux, we have to consider more 
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accurately the magnetic field structure in the GZK-sphere (see e.g Lemoine (2003) or Sigl 
et al. (1999)). 

An interesting consequence of the proton acceleration appears in considering the pp- 
collisions. Indeed, because n p ~ n e , the opacity of pp-collisions in a shell, namely t pp = 
n p (Tp P AR, is such that t pp = (upp/ar) r*. For this reason, we can easily write the radius r pp 
beyond which there is no more pp-collision: 



and r pp > r b for i] < which is quite possible according to the value of r)+ calculated in 
the previous section (this result slightly differs from which we found in Gialis & Pelletier 
(2003) because of a better estimate of r*). In this case, the proton acceleration starts in 
an opaque stage for pp-collisions: even if the energy limitation due to pp-collisions implies 
a cut-off energy around 1 GeV for the proton population, there is a possibility of a low 
energy emission of neutrinos between 5 GeV and 150 GeV which are produced via the ir + 
and the n~ decays (Gialis & Pelletier 2004). Moreover, the pp-process generates 7r°-mesons 
that decay in photons with an energy of 67 MeV in the co-moving frame. Because we have 
to consider the Klein-Nishina regime, the cross section, namely okn, is significantly lower 
than the Thomson cross section (<Jkn/&t — 1-7 x 10~ 2 ) so that the associated photospheric 
radius for these photons is about 0.13rv The photons produced at 67 MeV in a thin layer 
between 0.13 r* and r pp can cross over the medium without any electronic interaction; this 
is achieved if rj < ?7*/4. If one considers the pair creation process for 67 MeV photons on 
themselves, the cross section is such that <7 77 /<tt — 2.0 x 10~ 4 which leads to a transparency 
radius r 77 of about 10 _2 r + only. At last, these photons cannot interact with the thermal 
photons (< 0.5 keV) because the threshold energy is about 7.8 keV. 

Thus, we predict a possible signature of a such process which could be observed around 
20 GeV in the GRB spectrum when the baryon loading parameter is rj < 77*/ 4 which mainly 
occurs for a small number of shells (see Eq. (7)). At last, pp-collisions also generate a 
low energy emission of neutrinos because of the neutron decay which occurs after a proper 
neutron lifetime of ~ 13 hours. In fact, the neutron decay produces antineutrinos with 
energy about 0.7 MeV in the co-moving frame. This leads to antineutrinos of 0.2 GeV for 
the observer. 




(17) 



-10- 



3.2. Electron acceleration and energy limitation 



In this subsection, we analyse, the Fermi acceleration process concerning the electron 
population and we will use the same formalism than for the proton acceleration. Also, we 
assume that electrons and photons decouple at the photospheric radius r*. Thus, the electron 
acceleration starts at r + if r* > r b and at r b otherwise: we can define a radius corresponding 
to the beginning of the acceleration stage, namely r acc = max(r*,r&). In the co-moving 
frame, , the acceleration time for an electron with an energy e will be such that: 



t ~ 

^acc 



4 - dXiU ll0A300/ Wg) KlMeVJ U 



2 + OL 

3 



s. 



(18) 



We deduce that there are two main energy limitations in the electron acceleration beyond 
r acc . The first one is the synchrotron limitation which is: 

-5/4 



-syn 



~ 6.5 x 10 2 { — 



10 



300 



10 4 G 



(5q-2)/4 



MeV. 



(19) 



The second one is the expansion limitation: 



e exp ~ 1.2 x 10' 



\10J V300/ 



10 4 G 



1-0 



- MeV. 



(20) 



At rb, the synchrotron limitation is the strongest one but we can define a radius, namely r c , 
where these two limitations are equal and beyond which the main limitation is the expansion 
one: 

3/2 / £fo) N ~ 9 / 4 " 

,3ooy uo 4 ^ 



5.4 xlO- 5 f^) 9/4 f^- 
V10/ V30( 



(21) 



For a = 2, r c ~ 26 r b , but for a close to 1, this radius is beyond 10 5 r b . The electron energy 
increases and reaches the cut-off energy, e c , at the radius r c (see Fig. 3): this energy is very 
dependent on the magnetic field index a (B oc r~ a ). We have: 



e r ~ 6.5 x W 



10/ 



6a-3 



300 



(5.4 x 10~ 5 ) 



_ 5 \ 5a-2 



g(n) 
10 4 G 



-3 



MeV. (22) 



For a = 2, e c ~ 4.5 x 10 5 MeV and, for a close to 1, this energy reaches 1.2 x 10 7 MeV. 



For an electron with an energy e, the energy near which the synchrotron emission is a 
maximum can be expressed by: 

^ ,0 x 10 - gM) (£f (^) V, (23) 
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Fig. 2. — Left: Variation of r c /r^ as a function of the magnetic field index a for rj = 300. 
Right: Variation of the cut-off energy e c . 

and, according Eq. (19), is of the order of 0.3 keV at r b) in the co-moving frame, for accel- 
erated electrons around the synchrotron energy limit. 

Moreover, this energy reaches a maximum value at r c and we easily deduce that it is inde- 
pendent on a: indeed, it can be written in the co-moving frame: 



As we will see in the next section, this energy constitutes a cut-off energy in the high en- 
ergy emission spectrum, and turns out to be remarkably independent on the magnetic field. 
Contrary to the case of proton acceleration, the additional acceleration by the scattering off 
the hydromagnetic fronts is not operating because the transit time of the electrons accross 
any shell is too long compared to the synchrotron loss time for an energy above the estimate 
given by Eq.(19). 



In this section, we focus on the high energy emission due to both proton and electron 
population and we consider a magnetic field strength that decreases like r~ 2 . 




(24) 



4. Consequences on high energy emission 



- 12 - 



Log(e) 





(5a-2 )/4 

r 


j. 1- oc 




Synchrotron 


Expansion 


\v - 2/3 


limitation 


limitation 


Therm alisation 




1 <0C< 2 


y- e" 







1 1 — 

r b r* 



Fig. 3. — Diagram of energy limitation in the co-moving frame for the electrons. The position 
of the radius r c depends on the magnetic field index a and the electron energy at r c , e c varies 
from 10 5 to 10 7 MeV. 



4.1. Synchrotron emission by electrons 

We have seen in the previous section that electron acceleration is limited first by the 
synchrotron losses and, beyond r c , by the expansion losses. We can define the electron energy 
limitation, namely e&, by ef,(r) = e syn (r) up to r c and ef,(r) = e exp {r) beyond. 
According to Eqs. (23) and (24), there is a synchrotron emission over 3 energy decades 
between 0.3 and 200 keV in the co-moving frame, the highest energy corresponding to the 
cut-off energy hv c . Assuming an electron density, namely p(e,r), at the distance r and for 
the energy e such that 

p(e, r) oc e~ 2 exp (-e/ e 6 (r)) , (25) 
we deduce a local synchrotron spectrum i.e depending on the distance r which can be written: 

s e (z/,r)oc J p(e, r) P(e,r) 8{v - e 2 g(r)) de , (26) 

where P(e,r) oc e 2 g(r) 2 is the total radiated power for an electron and defining g(r) = 
e~ 2 i.e g(r) oc B(r) (see Eq. (23)). The integration of Eq. (26) easily gives: 

^r)oc 9 ( r )(^) I/2 exp(-^-|- y ). (27) 



V 
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Moreover, according to Eq. (23) we can define a local cut-off frequency, namely v+, by 
v*(r) = g(r) e b (r-) 2 so that we have v*(r) oc r 2 up to r c and v*{r) oc r -4 beyond. Thus, the 
integration of Eq. (27) leads to a spectrum: 



S e (u) oc v- 1/2 



r r " 3exp (-\/5y^) dr+ £ r ~ 3exp {-{^ r2 



- I dr 



{21 



where f*(r c ) corresponding to v c is the high cut-off frequency. 

Considering the case v < v c i.e around a few 0.3 keV in the co-moving frame, we show in 
appendix that Eq.(28) could be simplified and we obtain: 

Seiy) oc v- 1 . (29) 

Thus, electrons provide a gamma-ray emission which is constituted by a thermal component 
below a few 0.3 keV (see Sect. 2.2) and by a synchrotron component above. In the observer 
frame, the energy separating these two components could constitute the usual E peak (see 
e.g Band et al. (1993)) which is the energy for which the gamma-ray emission is the most 
important (i.e the peak energy of the burst). However most of the observations do not 
clearly exhibit such z/ 2 -spectrum in the low energy range. This E pea k is determined by the 
lowest value of the synchrotron cut off as displayed by the limitation energy diagram, which 
corresponds to an emission at r b . Moreover it depends on the baryon loading parameter rj 
because E pea k — rj h v syn {ri,) and can vary from a few 30 keV to about 300 keV in the observer 
frame. The energy spectrum of the gamma-ray emission will be such that 

5(1/) oc | V \ l T t^5T fe , • (30) 
[ v for E peak <hu<r]hu c 

With those estimates, there is a possible SSC generation under some conditions only: in 
fact, if r acc > r b , the electrons are quickly accelerated beyond 10 GeV and the SSC process, 
with keV photons, is in the Klein-Nishina regime (% h u syn 3> m e c 2 where 7 e is the average 
Lorentz factor of electrons). Because okn "C or, this leads to an associated photospheric 
radius smaller than r acc , so that there is no SSC generation. But, if r acc ~ rb, because of a 
strong synchrotron energy limitation for electrons, the SSC generation would give rise to a 
very low gamma emission in GeV range which would not be observable. 



At this stage, it is useful to remark that if we would have used a Bohm scaling, the 
electron energy cut off due to synchrotron emission would be much higher: 

isr^rfe)- <3i) 

This limit, that increases with distance contrarily to Kolmogorov one, would lead to a pro- 
hibitively high synchrotron spectrum. 
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4.2. Synchrotron emission by protons 

Let us first consider the acceleration of protons at internal shocks. The proton accel- 
eration via the usual Fermi acceleration is only limited by the expansion losses as we have 
seen in the Sec. 3.1. The previous calculation for electrons can be transposed to protons. 
For protons in the co-moving frame, the synchrotron emission is maximum at the energy 
(for Beer 2 ): 

^^.Oxxo-^Q- 2 (^V. ,32) 

In the observer frame, the resulting synchrotron spectrum for protons will be such that 
S p (u) oc u^ 1 . According to Eq. (14), this spectrum extends from 10~ 7 to 1 eV in the co- 
moving frame. These low energy photons can interact with the accelerated electrons via 
the inverse Compton scattering: in this case for which 7e h v <C m e c 2 with a high average 
Lorentz factor (% > 10 3 ) of the electrons, the cross section is a c ~ 7 2 a T . The associated 
photospheric radius, namely r p h, is such that (see Eq. (6)): r p h ~ 7 e rv 3> r b . During the 
internal shock phase, the inverse Compton effect thus produces photons with an energy am- 
plified by a factor 7 2 and which have an energy spectrum (oc e -1 / 2 ) extending from 0.1 to 
over 10 6 eV in the co-moving frame i.e from 10 eV to 0.1 GeV for the observer. 

Considering now UHECRs generated by the additionnal Fermi process with Lorentz 
factors, namely 7, in the range [10 8 , 10 9 ] in the co-moving frame (Gialis & Pelletier 2004): 
according to Eq. (16), the synchrotron energy spectrum will be such that S' u / iecr (z/) oc z/ -1 / 2 
and, following Eq. (32), UHECRs radiate synchrotron photons with an energy that scales 
like r~ 2 because of the magnetic field decreasing. Moreover, for a magnetic field decreasing 
like r -2 , the minimal value at r b must be > 10 5 G. According to the acceleration process, 
we have shown that a sizeable fraction of the UHECR component is not generated before a 
distance r cr of a few tens of r b typically between 50 r b and 100 r b and no longer after a few 
r cr . This UHECR generation radius r cr will determine the synchrotron emission range of 
UHECRs. In fact, because we have 

^,<r-10(^)(^ 2 (^)W, (33) 

the photons produced by UHECRs may extend from 1 MeV to 400 MeV in the co-moving 
frame. 

We have to examine, now, the consequences of a pair creation process between these 
hadronic photons and the electronic ones. We have previously seen that the energy of the 



- 15 - 



electronic photons reaches a maximum value of 200 keV at the distance r c that can be write 

r - 148 (iw) 3/4r6 ' (34) 

for B oc r~ 2 so that we have r cr < r c if r cr ~ 100 r b . According to the previous results (see 
Sect. 3.2) and the Eq. (33), the pair creation process happens if the product 



e p e e ~ 

syn syn 



™(^) 2 fe)" 4 (T^) 2 (^) 2( ^ )2 ' (35) 

is higher than 2 (m e c 2 ) 2 i.e ~ 5.2 x 10 5 (keV) 2 , where e e syn is the energy of electronic photons. 
The threshold Lorentz factor, j th , beyond which the hadronic photons undergo a pair creation 
can be defined by 

7 [ 2-3 x 10 8 (0) V4 for r cr <r<r c , 
|2.3xl0 8 (^)" 2 (^) 3 for r>r c . 

Thus, hadronic photons for which the Lorentz factor is between 10 8 and 7 S never undergo 
the pair creation process: so, the corresponding energy for hadronic photons will be in the 
range [^(7 = 10 8 ), e p syn (~f t h)} such that 

B(r b )\ ( r x 



eU7 = 10 8 ) ^ 1.0 l^g) (-) MeV, (37) 



z%n(lt h )^ 5.4 (^) (-) Me^, (38) 



and, 

'B(r b )\ 3/2 (^- 2 

with r cr ~ 100 r;, and for r cr < r < r c . Beyond r c , the cut-off energy for hadronic photons 

W ^ ^ 6 9/2 4 

(7*) ^2.5 f^Tl 9 2 f-) MeF. (39) 



10 5 G 



r, 



Therefore, the previous estimates allow to predict that an observer can detect syn- 
chrotron photons emitted by UHE-protons from a few 0.1 GeV to a few 10 GeV. We can 
remark that, if r cr > 225 r b then, e p syn {^ = 10 s ) < 200 keV, and this emission may not be 
observable because of the electronic synchrotron component. 



To end this review of all the possible absorption effects, we note that the interaction of 
these photons with themselves, which is kinematically possible, leads to a negligible opacity, 
namely, r 77 ~ (E syn /E) r* - (10" 5 - 10" 4 ) r*. 
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Let us estimate the corresponding global radiated energy that we will compare to the 
energy in the UHECR component. For an i? _2 -spectrum, assuming an uniform flux during 
At w , the energy E+ in the CR-component above 70 > 1 is 



„3 



E^ — TnipcQ r / p( 7 ) 7 d 7 At w , (40) 



7() 



where p( 7 ) = n*7 7 2 from 7 up to 7 max ~ 10 9 and the number of cosmic rays above 70, 
N+ = n*(rb) fir 2 c At w . We obtain 

^ = rm p c 2 iV, 7o log^^^) . (41) 

This can be simply compared to the GRB energy, E, since the total energy injected in 
protons, E p = -E*( 7 o = 1), is a sizeable fraction of E. We can write 

§~l- l0g7 ° , (42) 
E log 7 „ 



I max 



because N p ~ 70 iV*. Concerning the UHECRs with 70 = 10 8 , we reasonably find that 
E uhecr /E ~ 10" 1 . Let us come back now to the estimate of the radiated energy which could 
be observed. Beyond the generation radius, r cr , each proton synchrotron radiates a total 
energy 

which leads to: 

a_Y (0M\ 2 (_i_Y (i± 



«r(7)= 9.6X10'^) \^L) {^) \£) er g , (44) 

Thus, because dN uhecr /d'y — 10 s N uhecr 7~ 2 ~ N p ^~ 2 where N uhecr (7 > 10 8 ) and N p are 
respectively the total number of UHECRs and protons, the global emission which can be 
observed will have an energy: 



psyn 
uhecr 



/ iV^( 7 )7- 2 d 7 , ( 45 ) 
J 10 8 



and, with N p ~ E/(rjm p c 2 ), we obtain: 

For r cr = 50-100 r b) this ratio is between 10 -5 and 10 -4 . This leads to a number of photons of 
about 10" 3 -10~ 5 cm -2 for a GRB located at 1 Gpc and pointing towards observatory. But, if 
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we consider a magnetic field strength slightly higher of 10 5 G at r b , this number will increase 
around 10 -3 cm -2 . A few GeV-photons may be detectable by the GLAST instrument which 
will have an effective area of about 10 4 cm 2 . At last, this number of detected photons could 
increase up to about one hundred for a GRB located at 100 Mpc. 

This will constitute a very interesting signature of the UHECR generation in GRBs and 
will provide us some constraints on the internal shock model. The only possible competitive 
emission in that range could be produced by SSC process at the reverse shock and/or the 
early external shock under exceptional conditions (Granot & Guetta 2003). 

5. Conclusion 

The combined analysis of both electron and proton acceleration together with the ob- 
servation data leads to interesting conclusions in the frame of conservative assumption about 
the magnetic field and the scattering of particles off its irregularities. Indeed with a sub- 
equipartition magnetic field, that decreases like 1/r 2 , but concentrated in the shells invoked 
to account for the light curve, we have shown the following points. 

• The low gamma ray spectrum is satisfactorily explained by the synchrotron radiation 
of the electrons that are accelerated at the internal shocks with a spectrum displaying 
the expected peak emission. The lowest energy part of the spectrum would be likely 
explained by a thermal component. 

• The possibility of an UHECR generation under the previously stated conditions, as 
proposed in Gialis & Pelletier (2004), is confirmed and the estimated flux is in agree- 
ment with the expected one (Waxman 1995) to account for the Cosmic Ray spectrum 
around the ankle. A detailed estimate of the GRB contribution of the UHECRs to 
the Cosmic Ray spectrum, which is currently recorded by the Pierre Auger Observa- 
tory, has been proposed by Bahcall & Waxman (2001); it displays the expected excess 
around the GZK threshold. 

• The generation of cosmic rays in GRBs often starts before the fireball becomes trans- 
parent to pp-colisions. This gives rise to a gamma emission around 20 GeV due to 
7T°-decay. This emission is not contaminated by the SSC-emission of the electrons, be- 
cause that latter emission is in Klein-Nishina regime. We can also expect that neutrons 
are produced by these pp-collisions, they decay after some travel and thus generate a 
significant neutrino flux of 200 MeV (observer). 

• The most plausible signature of UHECR generation is not related to the p7-process, 
but more likely to their synchrotron emission. Indeed, the ratio of the corresponding 
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luminosities L pl /L syn ~ U s /U mag , where U s is the energy density of the soft photons 
and U mag the energy density of the magnetic field, turns out to be of order of unity. 
Since the energy of the neutrinos (~TeV for the observer) emitted through the P7- 
process is much higher than the energy of the gamma-photons (~GeV) emitted by 
the synchrotron process, the neutrinos are thousand times less numerous than the 
photons. Thanks to our conservative assumption about the magnetic field, a gamma- 
ray flux can be so emitted and could be observed. Indeed we found a natural range of 
gamma energies for which the fireball is transparent to pair creation process, because 
of the chronology of the emissions during the fireball expansion and of the decrease of 
the magnetic field in r~ 2 . This range typically extends from a few hundreds of MeV 
to 10 GeV. It turns out that the observation of GeV photons should be a signature of 
UHECRs with probably no confusion. This could be observed by GLAST. At 1 Gpc, 
a 10 5 G field would lead to a significant number of events which, of course, increases at 
shorter distances. This would be a very interesting signature of the UHECR generation 
in GRBs. 



A. Appendix 

We detail the exact determination of the synchrotron spectrum as defined by Eq. (28). 
Let us rewrite this expression: 



S e (v) CX V 



-1/2 



1 



— exp(-a— ) dr + / — exp(-o;— )dr 



r c 



(Al) 



where a = ^Jvjv c > y 'v syn {r b )/v c ~ 3.9 x 10 2 . 



First, a simple integration by parts leads to 

(a + 1) exp(— a) 

-ur L 

with A = rb/r c ~ 1/26 for B oc r~ 2 . 



1 r c I 

— exp(— a — ) dr = 



/ , 

J r r 3 r r^a 2 



a + A 
A 



exp(— a/X) 



(A2) 



For the second integral, because r c r^, we can assume r^ — > +oo and, next a first 
integration by parts, we can write for instance: 

I 7 , «P(-« ^ * = f a y ^4 - 4 J ? exp(-« -3) dr J (A3) 
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An iterative integration by parts gives, with some manipulations, the following result: 



1 r 2 exp(-a) 

— exp(-a — J dr = 



2 r 2 a 



E 

,n=0 



-l) n (n+l)! 



a' 



(A4) 



where the sum of the serie is convergent whatever a according to the generalized hypergeo- 
metric function (or Barnes's extended hypergeometric function). For a = a m i n ~ 3.9 x 1CT 2 , 
this sum is equal to 3.5 x 10~ 2 and for a = 1, the sum reaches 0.4. 
Thus, the synchrotron spectrum S e (u) is such that: 



S e (is) oc v 1 



a + 1 



a 



:-D"(»+i! 



n=0 



ft'' 



exp(— a) — a - exp(— a/A) 
A a 



(A5) 



For a m i n < a < 1, the expression between brackets is quasi constant so that S e (is) <x is l . 
Beyond a — 1, the same expression leads to a spectrum decreasing like v~ s with s £ [1.5, 2]. 
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